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INTRODUCTION 

Metal  matrix  composites  (MMCs),  because  of  their  favorable  elevated  temperature  properties,  wiii 
be  used  as  structure  materials  in  future  aircraft.  However,  at  elevated  temperatures,  these  composites 
lose  some  of  their  room  temperature  properties.  Moreover,  the  loss  becomes  much  worse  and  results 
in  material  deterioration  under  thermal  cycling. 

Thermal  fatigue  has  been  reported  to  occur  in  various  MMCs:  tungsten  fiber/copper  (1 ,2),  tungsten 
fiber/supperalloys  (3),  boron  fiber/aluminum  (4-6),  AljOj  (FP)  fiber/aluminum  (7),  FP  fiber/magnesium 
(8),  SiC  fiber/titanium  (9),  SiC  whisker/aluminum  (10),  and  graphite  fiber/aluminum  (11)  composites. 
Thermal  fatigue  is  caused  by  stresses  and  strains  due  to  repeating  constraints  of  free  thermal  expansion 
and  contraction.  The  constraints  can  be  grouped  into  two  categories:  external  and  internal.  The 
external  constraint  is  due  to  boundary  forces  applied  to  the  surface  of  the  component  which  is  being 
heated  or  cooled.  The  internal  constraint  is  produced  by  the  difference  in  coefficients  of  thermal 
expansion  (CTE)  of  the  constituent  phases  in  the  component  material:  CTE  mismatch,  and/or  non- 
uniform  temperature  distribution  in  the  component;  thermal  gradient.  CTE  mismatch-induced  thermal 
fatigue  is  a  characteristic  with  MMCs,  since  this  mismatch  for  most  of  MMCs  is  large. 

The  Coffin-Manson  relationship  for  low  cycle  fatigue  (12,13)  is  reported  to  be  applicable  for  thermal 
fatigue  (13,14).  Therefore,  the  thermal  strains,  elastic  and  plastic,  generated  by  the  CTE  mismatch,  are 
essential  in  the  analysis  of  thermal  fatigue  of  MMCs.  Thermal  stress  and  strain  have  been  anaiyzed 
extensively  for  continuous  fiber  (1 5-31 )  and  whisker  reinforced  (32, 33)  MMCs.  However,  such  analysis 
has  been  very  limited  for  particulate  MMC  (34,  35). 

This  report  presents  the  analysis  of  thermal  stresses  and  strains,  caused  by  CTE  mismatch  between 
the  spherical  reinforcement  particle  and  the  matrix  in  an  MMC,  and  the  estimation  of  the  thermal  fatigue 
iife.  In  the  analysis,  the  following  assumptions  are  used. 
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1 .  The  reinforcement  particle  is  a  sphere  in  an  Infinite  matrix  of  an  MMC. 


2.  The  matrix  has  elastic-perfectly  plastic  behavior. 


3.  The  stress-strain  behavior  is  independent  of  strain  rate  and  stress  orientation. 


4.  The  temperature  in  the  MMC  is  uniform  at  all  time.  Therefore,  there  Is  no  thermal  gradient  in 
the  MMC. 


STRESS  AND  STRAIN  ANALYSIS 

The  thermal  stresses  and  strains,  induced  by  CTE  mismatch,  will  be  analyzed  under  purely  elastic 
conditions  first  and  then  under  plastic  conditions. 

A.  Thermoelastic  Stresses  and  Strains 

The  following  elastic  model  was  used  in  the  artalysis  of  thermoelastic  stresses  and  strains. 

1 .  At  high  temperatures,  the  radii  of  the  matrix  hole  and  the  reinforcement  particle  are  identical  R, 
(Figure  1). 

2.  Upon  cooling,  the  matrix  hole  tries  to  contract  to  a,  and  the  reinforcement  particle,  with  a  lower 
CTE.  tries  to  contract  only  to  r^,  (Figure  1). 

3.  The  final  or  effective  radius  of  the  reinforcement  particle  under  the  constraint  of  the  matrix  is  r,. 
(Figure  1). 
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The  radii  arKJ  their  relationship  can  be  denoted  by  the  following  equations. 


a  =  R(1  -  «„AT) 

(1) 

ro  =  R(1  -  OpAT)  =  a(1  +«) 

(2) 

r,  =  a(1  +  0S) 

(3) 

where 

:  CTE  of  matrix 

«p  :  CTE  of  reinforcement  particle 

s  misfit  factor 

0  constraint  factor 

AT  ;  range  of  temperature  fall 

Aa  : 

From  Equations  (1)  and  (2),  the  misfit  factor  s  can  be  described  as  a  function  of  difference  in  CTEs, 
range  of  temperature  fall  AT,  and  CTE  of  matrix  Equation  (4). 

« -A./[  {1/AT)-a„]  (4) 


Throughout  this  analysis,  spherical  coordinates  are  used  The  origin  is  located  at  the  center  of  the 
spherical  reinforcement  particle.  Because  of  symmetry,  the  tangential  displacement  as  well  as  the  shear 
stress  and  strain  are  all  zero,  and  the  radial  displacement  u  is  a  function  of  radial  distance  r.  There  are 
three  non-zero  stress  components,  a  radial  stress  a,  and  two  tangential  stresses  a,  and  a^.  These 
stresses  must  satisfy  the  equilibrium  condition  in  the  radial  direction,  Equation  (5). 
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do,/dr  +  (2/r)  (0,-0,)  -  0  (5) 

Furthermore, 

o,  -  o,  and  <•  *  <4  (6) 

In  a  cooled  body,  the  total  strain  Is  made  up  of  two  components  in  each  of  the  radial  and  tangential 
directions.  Equations  (7)  and  (8).  One  component  is  a  uniform  contraction  proportionai  to  the  range 
of  temperature  fall.  Since  this  contraction  is  equal  in  di  directions,  only  normal  strains  and  no  shear 
strains  are  found.  Thus,  the  normal  contraction  in  any  direction  is  •  «AT.  The  other  component 
comprises  the  strains  required  to  maintain  the  continuity  of  the  body.  These  strains  are  related  to  the 
stresses  by  means  of  Hooke's  law  of  isothermal  elasticity. 


€,  -  (1/E)  (Or-avo,)  -  CAT  (7) 


-  (1/E)  •  (  -  VO,  ♦  (  1  -  V  )  o,  1  -  oAT  (8) 

Furthermore,  the  radial  and  tangential  strain  components  are  defined  as 

e,  -  du/dr.  c,  -  u/r  (8) 


where 

r :  radial  distance 

a„  t, :  radial  stress  8:xl  strain 

o*.  04.  :  («,  <4  :  tangential  stresses  arxJ  strains 

u  ;  radial  displacement 
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V  Poisson’s  ratio 

E  elastic  modulus 

o  :  CTE 

From  Equation  (5)  and  Equations  (7)  -  (9),  the  equilibrium  equation  can  be  expressed  with  radial 
displacement  u  and  distance  r. 


d*u/dr*  +  (2Jr)  •  (du/dr)  -  2u/r*  -  0 


(10) 


The  general  solution  of  Equation  (10)  is: 

u  -  C^r  +  Cj/r* 


(11) 


The  boundary  condition  for  the  matrix  is 

(uj,.,,-r,  -a-ap&.  (uJ,..-0  (12) 


where  u„  is  the  radial  displacement  in  the  matrix. 

Then, 

C,  -  0,  Cj  -  a®p5  (Wp«) 
(uj,,„  -a’p»  (1  *  Pft)/r* 


(13) 


The  bourxJary  condition  for  the  spherical  reinforcement  particle  is 

(Up),-r,  -  ri  -  r,  -a(P  -  1)5.  (Up),.,-  0 
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where  Up  is  the  radial  displacement  in  the  spherical  reinforcement  particle. 
Then, 


C,  -  O  -  1)5 /(I  P5).  Cj-O 

(Up),<,,  -  (P  -  1)»r/(1  ♦  P5) 


Substituting  this  value  of  Up  into  Equation  (9)  gives  the  elastic  strains  within  the  reinforcement  particle. 

(er)p  -  (ee)p  -  (P  -  1)8 /(I  ^  P6)  (15) 


From  Equations  (6)  -  (8)  and  Equation  (15),  the  elastic  stress  components  within  the  reinforcement 
particle  are  determined  to  be 

(°r)p  "  (®e)p  “  (®*)p  ^5) 

-  [  Ep/(1  -  2vp)  1  ■  [  (p  -  1)5  /  (1  *  P5)  .  apAT  1  ^  ' 

This  equation  indicates  that  the  spherical  reinforcement  particle  is  in  a  state  of  hydrostatic  stress,  which 
is  larger  with  larger  elastic  modulus,  GTE,  and  temperature  range.  Consequently,  the  reinforcement 
particle  is  in  an  elastic  state  and  it  is  not  plastically  deformed. 

Substituting  the  displacement  value  in  the  matrix.  Equation  (13),  into  Equation  (9)  furnishes  the 
elastic  strain  components  in  the  matrix  as  follows. 


(€,)„  -  -  2  (€,)„  -  -  2P5(1  *  P5)  -(a/r)® 

-  -  Hr,  I  rf 


(17) 
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where 

X  -  2pS(1  +  ^ft)‘ 

The  feature  of  elastic  strain  component  variation  with  radial  distance  In  the  matrix  is  shown  in  a  pIcM  of 
elastic  strain  component,  normalized  with  X,  vs.  radial  distance  r,  Figure  2.  The  magnitude  of  each 
elastic  strain  component  is  largest  at  the  reinforcement  partide/matrix  interface,  r  -  r„  arxJ  decreases 
with  increasing  distance  from  the  interface. 

From  Equations  (7),  (8)  and  (17),  elastic  stress  components  in  the  matrix  are  fourxf  to  be 


(a,)„  -  -  [  2E„p«(Upa)  /(I  («/r)>  *  (E„a„AT)  /  (1  -  2v„) 

-  -  2Q„X(r,/r)®  .  3B„«„AT 


(18) 


(c,)^  -  I  E„pa(1  >  p«)/(1  *  V  J  1  •  (a/r)»  .  (E„«„AT)  /  (1  -  2v„) 
-  Q„A(r,/r)®  ♦  3B„«*AT 


(19) 


where 

G„  :  shear  modulus  of  matrix 
B„  :  bulk  modulus  of  matrix 

The  features  of  the  elastic  stress  component  variation  with  radlai  distance  in  the  matrix  is  shown  in  a 
plot  of  [  (<r,)„  -  3B„  a„,  AT  ]  /G„X  VS.  r  and  [  -  3B„  m„  AT  ]  /G„X  vs.  r.  Figure  3.  The  magnitude 

of  each  elastic  stress  component  Is  also  largest  at  the  reinforcement  partide/matrix  intertace,  r  >  r„  arxf 
decreases  with  increasing  distance  from  the  interface. 
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B.  Plastic  Deformation 

As  was  pointed  out,  the  spherical  reinforcement  particle  rerruiins  in  an  elastic  state  and  the  stresses 
are  greatest  in  the  matrix  adjacent  to  the  reinforcement  particle.  For  the  analysis  of  plastic  deformation, 
the  following  model  is  taken. 

1 .  A  thick  hollow  sphere  is  deformed,  elastically  and  plastically,  under  uniformly  distributed  internal 
pressure.  Figure  4. 

2.  The  internal  pressure  is  prcxluced  by  a  misTitting  spherical  reinforcement  particle. 

The  internal  pressure  p  is  identified  as  the  radial  stress  component  at  the  reinforcement  particle/matrix 
interface,  r  =  r,.  Thus  from  Equation  (18), 

P  -  (o,)„  -  -  [{2E„pa(1-P»)  J/(1  >  vjl  .(a/r,)»ME„,«„AT)/(1-2vJ  ^go) 

-  -  2Q^X  .  3B,«„AT 

According  to  the  von  Mises  yield  criterion  (36),  plastic  deformation  occurs  if  an  equivalent  stress 
((7,)„,  defined  by  the  foliowing  formula,  reaches  the  yield  stress  of  the  matrix 

(o,)„  -  (l/v/2)  .  I  (  (a,)„  -  (o,)„  )2  .  {(a,)„  -  (o^)„  }  *  .  ^21) 

Since 

(Oe)m  -  (sim.  (o.)m  ’  (o*)m  "  (®r)m  (22) 
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From  Equation  (19)  with  r  =  r,,  Equation  (20),  and  the  above  equation,  plastic  deformation  starts  at  the 
reinforcement  particle/matrix  interface,  if  the  following  condition  is  satisfied. 

(o.)„  -  -  (3/2)  Ip,  -  (E„a„AT)/(1-2v  J1  -  (a,)„  (23) 

where  p,  is  the  critical  internal  pressure  for  plastic  deformation  (or  yielding)  of  the  matrix. 

From  Equation  (23), 

p,  -  -  I  (2/3)(Oy)„  -  (E„a„AT)  /  (1  -2v  J  )  (24) 


This  equation  indicates  that: 

1 .  The  critical  internal  pressure  for  plastic  deformation  is  less  than  2/3  of  the  yield  strength  of  a 
given  matrix  material. 

2.  Its  magnitude  decreases  linearly  with  increasing  range  of  temperature  fall  and  decreasing  yield 
strength. 

3.  Its  magnitude  is  reduced  further  for  any  other  matrix  material  with  larger  elastic  modulus  and 
CTE. 

With  increasing  internal  pressure  beyond  the  critical  value  p^,  a  plastic  zone  is  formed  in  the 
immediate  vicinity  of  the  reinforcement  particle/matrix  interface  and  extended  to  a  certain  radius  r,.  This 
radius  r,  separates  the  inner  plastic  zone  from  the  outer  elastic  zone.  The  equilibrium  equation  for  the 
plastic  zone  is  obtained  by  substituting  the  yield  condition 
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(®«)m  “  (®y)m  “  ~  (®f)m 


into  Equation  (5).  The  equilibrium  equation  is 


d(o,)  Jdr  -  2(0,)  Jr  -  0 


(25) 


Integrating  Equation  (25). 


(o,)„  -  2(0, )„  •  In  r  +  C 


(26) 


Using  the  boundary  condition  (o,)„  =  -  p  at  r  =  r,,  the  constant  C  becomes 


C  -  -  p  -  2(o,)„  •  In  r. 


(27) 


From  Equations  (25)  -  (27).  the  stress  components  in  the  plastic  zone,  (r,  <  r  <  r^).  are 


(Or)m  “  2(o,)„  Nr/r,)  -  p 


(28) 


-  (o,)„  [2  Nr/r,)  *  1]  -  p 


(29) 


The  stress  components  in  the  elastic  zone,  (r^  <  r),  can  he  found  by  substituting  p^  for  p  and  r^  for  r, 
in  Equations  (28)  and  (29)  as  follows. 

(o,)„  -  (o,)„  I  (2/3)  -  2  •  ln(r/r,)  J  *  (E„a„AT)/(1  -  2v„)  (30) 
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-  (Oy),n  I  (3/5)  *  2  •  ln(r/r,)  I  -  (E„a„AT)/(1  -2v„)  (31) 

At  the  boundary,  which  separates  the  plastic  and  elastic  zones,  the  radial  stress  components  are 
identical,  I.e.,  Equation  (28)  =  Equation  (30)  at  r  =  r^.  From  this  relationship,  the  plastic  zone  radius 
r^  is  determined  as 

r,  -  r,  •  exp(  (5/6)  .(1/(2  (o,)  J  }  •  (p  -  (E„a„AT)/(1  -2v„)  }  ]  (32) 

By  substituting  of  Equation  (30)  and  of  Equation  (31)  into  the  equations  of  stress-strain 
reiation.  Equations  (7)  and  (8),  the  radial  and  tangential  strain  components,  («,)„  and  (eg)„.  in  the 
elastic  zone  of  the  matrix  is  found  as 

(Cr)™  ■  (<'y)m/Em)  *  I  (2/3)  *  (1-5v  J  .  (1-2v„)  •  ln(r/r,)*  ]  ,33. 

-2a„AT 

(€e)m  -  •  I  (1/3)  •  (5-7v„)  .  (1  -2v  J  •  ln(r/r,)2  J  (34^ 

-2a„AT 

In  the  plastic  zone,  the  stress-strain  relations  are  described  as 

(€,)„  -  du/dr  -  (1/E„)  •  [  (o,)„-2v  Ja,)„  J-a„AT.(€,)„P  (35) 

(ee)„  -  u/r  -  (1/EJ  •  I  -vJo,)„.(1  -V  J  •  (o.)„  ]-a„AT.(€«),„P  (36) 
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where  («,)„'’  and  («»)„’’  are  plastic  strain  components  in  radial  and  tangential  directions,  respectively. 
From  the  incompressibility  condition  for  plastic  strains.  («,)„”  +  2  («,)„”  =  0,  and  Equations  (35)  and 
(36),  the  following  differential  equation  is  derived. 

du/dr  +  2u/r  -  (1/B„)  [  2(o^)„  •  ln(r/r,)  .  (  2(Oy)„  -  3p  )  /3  ]-3a„AT  (37) 


where  =  E,„/3(1  -  2v  J  is  the  bulk  modulus  of  the  matrix.  Its  general  solution  is 


u  -  2(oy)„  •  r  ■  ln(r/r,)/3B„-pr/3B„-a„  •  r  •  AT^C/r^ 


(38) 


The  constant  C  is  calculated  from  the  displacement  at  the  plastic-elastic  zone  boundary  u,  given  by 
Equation  (34),  for  r  =  r<j  as  follows; 


C  -  (Oy)m 


(39) 


From  Equations  (38)  and  (39),  the  displacement  in  the  plastic  zone  (r,  <  r  <  r^)  u  can  be  written  as 


u  -  2(a,)„  •  r  •  ln(r/r,)/3B„  -  pr/3B„  *  (  (0y)„v  )  •  (r.^’/r*) 


-  a „  •  r  •  AT 


(40) 


Knowing  the  displacement  u.  the  strain  components  in  the  plastic  zone  (r,  <  r  <  r^)  can  be  determined 
as  follows. 


(€,)„-du/dr  -  {  2(oy)  J3B„  )  •  ( ln(r/r,)  .  1  }-p/3B^ 
-(2(o,),  vJE„}  (r,/r)3-a„AT 


(41) 


(€e),-u/r  -  {  2(0,)  J3B„  }  •  (ln(r/r,)- p/3B„ 
-{(a,)„v„/E„)  (r,/r)3-a„AT 


(42) 
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Employing  the  stress-strain  relations,  given  by  Equations  (35)  and  (36),  the  yield  condition  {o^)^  = 

-  ((7,)„„  and  the  incompressibility  condition  (£,)„"  the  plastic  strain  components  are 

determined  as 

(£,)„»  -  -2(€«)„p  -  {  2(a,)  }  -  I  1  -  {  1  .  (r,/r)3  }  ]  (43) 

The  elastic  strain  components,  («,)„*  and  (c,)^*  are  obtained  by  subtracting  the  plastic  strain 
components.  Equation  (43),  from  the  total  strains.  Equations  (41)  and  (42). 

(er)m*  -  («r)m  '  (er)m‘’ 

-  (1/3BJ  •  {  2(o^)„  •  ln(r/r,)  -  p  }  -  2(o^)„v  JE„  (44) 


(ee)m*  -  (€«)„,  -  (es)^'’ 

-  (1/3BJ  •  <  2(a,)„  •  ln(r/r,)  -  p  )  -  {  (o,)„/EJ  •  (1  -  v  J  }  (45) 


THERMAL  FATIGUE  LIFE 


On  the  basis  of  Manson's  (13)  and  Garmong’s  (14)  reports  on  the  applicability  of  the  low  cycle 
fatigue  damage  model  to  the  thermal  fatigue,  the  Coffin-Manson  relationship  (12,  13)  can  be  utilized  for 
the  estimation  of  thermal  fatigue  life. 

(€)P  -  M.N,"  (46) 


where 

(f)® :  plastic  strain  per  thermal  cycle 

M,  z  ;  material  constants 
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N,  :  number  of  thermal  cycles  at  failure. 

With  total  strain,  Manson’s  simplified  equation  (37)  may  be  employed. 

(c)  -  3.5  •  (oJE)  ■  •  (N,)-°« 


where 

(«)  :  total  strain  per  thermal  cycle  =  («)*  +  («)’’ 
ultimate  strength 

(c,)  :  true  fracture  strain  in  a  monotonic  tensile  test 


(47) 
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CONCLUSIONS 

1 .  In  a  MMC  with  a  spherical  reinforcement  particle,  the  CTE  mismatch-induced  thermoelastic  stresses 
and  strains  are  largest  at  the  reinforcement  particle/matrix  interface,  and  they  decrease  with  distance 
from  the  interface. 

2.  The  spherical  reinforcement  particle  is  in  a  state  of  hydrostatic  stress,  which  is  larger  with  larger 
elastic  modulus,  CTE,  and  temperature  range. 

3.  Plastic  deformation  starts  in  the  matrix  adjacent  to  the  reinforcement  particle.  Consequently,  the 
interface  is  a  potential  site  for  crack  initiation  under  thermal  cycling. 

4.  The  critical  internal  pressure  for  plastic  deformation  is  less  than  2/3  of  the  yield  stress  of  a  given 
matrix  material,  and  it  decreases  with  increasing  range  of  temperature  fall. 

5.  With  the  analytically  determined  strain,  the  thermal  fatigue  life  can  be  estimated  by  employing  the 
Coffin-Manson  relationship. 
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Figure  3.  Elastic  Stress  Cornponent  Variation  with  Radial  Distance. 
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.  Thick  Hollow  Sphere  Under  Internal  Pressure. 
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